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ABSTRACT 

We investigate the implications for some nonstandard cosmological models 
using data from the first three years of the Supernova Legacy Survey (SNLS3), 
assuming a spatially fiat universe. A comparison between the constraints from 
the SNLS3 and those from other SN la samples, such as the ESSENCE, Union2, 
SDSS-11 and Constitution samples, is given and the effects of different light-curve 
fitters are considered. We find that SN la with SALT2 or SALT or SIFTO can give 
consistent results and the tensions between different data sets and different light- 
curve fitters are obvious for fewer-free-parameters models. At the same time, 
we also study the constraints from the SNLS3 along with data from the cosmic 
microwave background and the baryonic acoustic oscillations (CMB/BAO), and 
the latest Hubble parameter versus redshift (H(z)). Using model selection crite- 
ria such as x^/dof, GoF, AIC and BIC, we find that, among all the cosmological 
models considered here (ACDM, constant w, varying w, DGP, modified poly- 
tropic Cardassian, and the generalized Chaplygin gas), the fiat DGP is favored 
by the SNLS3 alone. However, when additional CMB/BAO or H(z) constraints 
are included, this is no longer the case, and the fiat ACDM becomes preferred. 

Subject headings: cosmology:observations - supernova: general 
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INTRODUCTION 



Type la supernovae (SNe la) have altered the focus of cosmology dram a.tically since 



they f irst indicated that the ex pansion of the Universe is currently accelerating ( iRiess et al. 



19981 : iPerlmutter et al. Ill999l ). Usually, the mysterious cause of the cosmic acceleration is 



attributed either to the existence of an exotic energy component, called dark energy, or a 
modification of the standard theory of gravity. At present, SNe la are still the most direct 
probe of the history of cosmic e xpansion and the propert ies of dark energy, and many SN 



la samples, s uch as ESSENCE fWood-Vasev et al. 1120071). SDSS-II (Kessler et al. 1 120091). 



Constitution (IHicken et al. II2009I ) and Union2 (lAmanullah et al. Il2010l ). have been released. 



The Supernova Legacy Survey (SNLS) is an ongoin g five year project th at aims to 



probe the expansion history of the universe using SNe la ( Sullivan et al. II2005I ). The goal 
of this survey is to measure the time-averaged equation of state (EOS) of dark energy {w) 
to 5% (statistical uncertainties only) in combination with other cosmological probes and to 
10% including the systematic uncertainties. Analys is of the first year o f data from SNLS, 
containing 71 high-redshift SNe la, was presented in lAstier et al. I (120061 ). 



Recently, data f rom the first three years SNLS data set were released. iGuy et al. 
(I2OIO. hereafter GlOl ) presented the photometric properties and relative distance moduli of 
252 SNe la from this sample. Based on selection criteria, they removed 21 SNe and obtained 
Vim = 0.211 ± 0.034(stat) ± 0.069(sys) for a flat ACDM using the SNLS3 data alone. Using 



a slightly different set of selection criteria, IConley et al. I ( 12011. hereafter Cll ) selected 242 



SNe la from SNLS, and combined these with 123 low-^, 93 intermediate-^ (SDSS-11) and 14 
high-z SNe from the Hubble Space Telescope to form a high-quality unified sample containing 
472 SNe la. We shall refer to this sample as SNLS3. For the SNLS3 data alone, Cll finds that 



w 



-0.91_o;2o(stat)_o;i4(sys), considering a fiat Universe with constant w. The constraints 
on dark energy w hen SNLS3 is combin ed with other probes (such as the b aryonic acoustic 



m 



Sullivan et al. I ( 



os cillation (BAG) ( Percival et al. 1120091) or WMAP7 (IKomatsu et al. II2OIII )) were presented 



2011. hereafter Sill ). The results for a fiat constant w universe are = 
-1.0 611nnfi«i and those for a flat varying w uni verse characterized by 
the CPL parameterization ( IChevallier &: Polarski II2OOII : iLinder Il2003[ ) (w(a) = wq + Wa{l — 



0.269 ±0.015 andw 



a)) are wq = -0.905 ± 0.196 and Wa = -0.9841};°^^. 

Dark energy is not the only possible explanation of the present cosmic acceleration, 
and many other models have been considered such as DGP, Cardassian, and the generalized 
Chaplygin gas. In this paper we investigate co nstraints on t h ese no nstan dard cosmological 



model s. Similar analyses we re carried out bv iDavis et al. I ( 120071 ) and 



(]2009f ) using the ESSENCE (IWood-Vasev et al. 1120071 ) and SDSS-II (IKessler et al. 1 120091 ) 



SoUerman et al. 



SN la data sets, respectively. In order to give a comparison between the results from SNLS3 
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with those from other SN la data sets, we also compare the constraints obtained using SNLS3 
to those using other popular SNe la samples such as the ESSENCE, Constitution, Union2 
and SDSS-II samples, and additionally analyze the effects of different light-curve fitters on 
the results. In addition to the SNe la data, constraints from other probes, such as the cosmic 
microwave background (CMB) and baryonic acoustic oscillations (BAO) measurements of 
the acoustic scale (CMB/BAO) and the Hubble parameter vs. redshift (H(z)), are also 
considered. 

We organize our paper as follows. Section 2 presents the data sets and the statistical 
analysis method. Section 3 describes briefly the models considered and the observational 
constraints on them. Model testing results using model selection statistics are shown in 
section 4. Finally, Section 5 gives a summary of our results. 



2. DATA SETS AND ANALYSIS METHODS 



In this Section, we describe the data sets and analysis techniques. 



2.1. Type la supernovae 



The SNLS3 SN la sample has been discussed in detail in Cll. This sample consists 
of 472 SNe la (123 \ow-z, 93 SDSS, 242 SNLS, and 14 H ubble Space Telescope) and is 
compiled with the combined SiFTO JConlev et al. I boosh and SALT2 JGuv et al. I l2007h 
light curve fitters. As discuss ed in Sll, the SNLS 3 data set has several advantages over the 
first year sample presented in lAstier et al. I (120061 ). First of all, the sample size has increased 



by a factor of three. Furthermore, the sources of potential astrophysical systematics are 
examined by dividing the SN la sa mple according to the properties of either the SN or its 



environment (jSuUivan et al. Il2010l ). Finally, an ameliorated photometric calibration of the 



light cur ves and a more consi stent understanding of the experimental characteristics are 



allowed (IRegnault et al. II2009I ). 



For the cosmological analysis, one can minimize the following 



X 



E 

SN 



stat 



int 



where is the rest-frame S-band magnitude of an SN, 



m 



mod 
B 



is the magnitude of the SN 



predicted by the cosmological model, a stat includes the uncertainties both in and m^"'^ 
and aint parameterizes the intrinsic dispersion of each SN sample. The values for a^t used 
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in our analysis can be found in Table 4 in Cll. The mo del- dependent magnitude is given by 
m^°'^ = 5 logio ^^L(^hci, 2;cmb; p) - a{s - 1) + /3C + Mb, (2) 

where p stands for the model parameters set, Vl is the Hubble-constant free luminosity 
distance, Zhei and Zcmb are the heliocentric and CMB frame redshifts of the SN, and a and 
P characterize the stretch and color (s and C)-luminosity relationships. Ais represents 
some combination of the absolute magnitude of a fiducial SN la and the Hubble constant. 
Here, a, /3 and Mb are "nuisance parameters". We marginalize over A^s by following the 
method described in Cll (Appendix C). In this marginalization method the environmental 
dependence of SN properties, which allows A^^ to be split by host-galaxy stellar mass at 
1O^°M0 (this correction is not employed for all other SN la samples we are considering), is 
taken into account. Defining the vector of residuals between the model magnitudes and the 
observed magnitudes Am = — m^°'^, we obtain the observational constraints on model 
parameters by minimizing 

= AA'^ • C-i • Am, (3) 

where is the inverse covariance matrix [^. The total covariance matrix is detailed in 
section 3.1 in Cll. The covariances which take both statistical and systematic uncertainties 
into consideration are used in our investigation. We allow a and /3 to vary with the cos- 
mological parameters to avoid biasing our results using a grid minimization routine and 
then use their best fit results to obtain the allowed regions of model parameters. 

It is well-known that to obtain the precise distance of the standard candle, SNe la, is 
crucial for modern cosmology. Since the intrinsic luminosity of SNe la is correlated with 
the sh ape of its optic al light curves, one can determine it by proposing a method to relate 



them. iPhillips I (119931 ) firstly introduced a method by finding a correction between the SNe la 
intrinsic luminosity and the parameter Amis (5), where Amis (5) is the amount of a SNe la 
B-band declination during the first fifteen days after maximum light. The MLCS is a different 
m ethod using the mu lticolor light curve shapes to estimate the luminosity distance proposed 



m 



RiessetaL 



Jha 



( 1995). Th is approach was extended to include U-band measurements, named 



(I2OO2I ). The MLCS2k2 has been applied for distance estimate in several 



MLCS2k2, by 

popular SN la samples, such as ESSENCE, SDSS-II, and Constitution. 



Guy et al. I (120051 ) proposed another method, SALT, to construct the SNe la luminosity 
distance estimate by parameterizing the light curve with a parameter set, i.e., a luminos- 
ity parameter, a decline rate parameter and a single color parameter. It offers a practical 



^https:/ /tspace. library, utoronto.ca/handle/1807/26549 
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advantage w hich makes i t easi ly apphcable to high-redshift SNe la. By including spectro- 
scopic data, iGuy et al. I ( 120071 ) improved this SALT to SALT2. SALT2 has been applied 



to calculate the distance modulus for several popular SN la samples, such as Constitution, 
Union2, SDSS-II, and SNLS3. It is notable that the SDSS-II corrected the SALT2 results for 
selection biases by using a Monte Carlo simulation (see section 5.2 and section 6 of Kessler 
et al. 2009) and the SNLS team made a few technical modifications in the training proce- 
dure, such as higher resolution for the components and the color variation law, and a new 
regularisation scheme, to S ALT2 for the SNLS3 (detailed in section 4.3 and Appendix A of 



Guy et al. 2010). Recently, IConley et al. I ( 120081 ) proposed another new empirical way, i.e., 
SiFTO, which is similar to SALT2 using the luminosity magnitude, light curve shape and 
color to obtain the distance. This method has been applied in the analysis of SNLS3. 



P revious works have shown that different SN la data sets give different results (IRydbeck et al 



20071 : iBueno Sancheza et al. I l2009l ) and different light-curve fitter s may also lead to dif- 



ferent conclusions although the same SN la sample is considered ( ISoUerman et al. 1 12009 



Bengochea I l201ll ). Therefore, in the present paper, besides the SNLS3 (analyzed with 
SALT2 and SiFTO light-curve fitters), we also consider some other popular SN compilations: 
ESSENCE (SALT and MLCS2k2), Constitution (SALT2 and MLCS2k2 with Ry = 1.7), 
Union2 (SALT2) and SDSS-II (SALT2 and MLCS2k2). However, it should be mentioned 
that these SN la samples are not independent, since many of them are drawn from the same 
sources. That is , dif ferent SN la com pilations m ay contain many of the same SNe. For 
instance, the nearby (jjha et al. 1 120071 ) and SNLS ( lAstier et al. I l2006l ) SN la are, in part, 
included as subsets in all above mentioned samples. In addition, except for the Union2 and 
SNLS3, there is no easy way to include systematic uncertainties for other referred SN sam- 
ples. In our analysis only for the SNLS3 sample are the systematic uncertainties considered. 



2.2. The Cosmic Microwave Background and Baryon Acoustic Oscillations 

The CMB and BAO constraints used in this analysis are based on the angular scales 
measured at the CMB decoupling epoch at z^, = 1090 and imprinted in the clustering of 
luminous red galaxies (BAO) at z = [0.2,0.35]. Its acoustic scale is given by 

dA{z^){l + z^) 

Ia = TT , (4) 

which represents the angular scale of sound horizon at decoupling, where = dL/{l + z)"^ 
is the proper (not comoving) angular diameter distance, di = (1 + -2) ^ is the luminosity 
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distance, and rs(z*) is the comoving sound horizon at recombination, 



Cs[Z) 

H(z) 



dz 



(5) 



which depends on the speed of sound, c^, in the early universe. The BAO scale is given by 
rs{zd)/Dv, where the rs{zd) is the comoving sound horizon at drag epoch {zd ~ 1020), and 
the so-called dilation scale, Dy, is a combination from angular diameter distance and radial 
distance according to 



;i + zfd' 



cz 



1/3 



(6) 



^By combining the BAO measurements of rs{zd) / Dy{z) at two redshifts (jPercival et al. 

2007h . rs{zd)/Dv{Q.2) = 0.1980 ± 0.0058 and r,(zd)/ D y (0.35) = 0.1094 ± 0.0033, with the 
CMB measurement of Ia given in lKomatsu et al. I ( l2009l ). Ia = 302.10±0.86, and considering 
the ratio of the sound horizon at the two epochs, rs{zd)/rs{z^,) = 1.044 ± 0.019, the final 
constraints we use for the cosmology analysis are obtained: 

dA{z^){l + z^) 



/o. 



20 



0.35 



Dv{z = 0.2) 
dA{z*){l + z*) 



19.04 ±0.58 



10.52 ±0.32, 



(7) 



Dv{z = 0.35) 

which do not depend on the comoving sound horizon scale at recombination. However, these 
two ratios are not independent with a correlation coefficient of 0.39. This correlation is 
considered in our cosmological fits. Moreover, the constraints from these two ratios (labeled 
CMB/BAO here) are expected to be a good approximation for all models tested in this paper 
since the redshift difference between the decoupling and the drag epoch is relatively small, 
and the sound horizon at these two epoch is mostly dominated by t he fractional difference 
between the number of photons and baryons (jSoUerman et al. 1 120091 ). 



In addition, we do not use the CMB "shift parameters" {Ia, R, 2;*) derived from WMAP7 (IKomatsu et a 



20111 ) in our analysis, because the CMB distance prior is applicable only when the model 
in question is based on the standard Friedmann-Lemai tre- Robertson- Walker universe with 
matter, radiation, dark energy and spatial curvature (IKomatsu et al. II2OIII ). and we are 



testing some nonstandard models, e.g. DGP and modified polytropic Cardassian. 



2.3. Hubble parameter versus redshift data 



The Hubble parameter versus redshift data to be used in this paper consists of three 
subsamples. The first one includes 9 data points which are derived from differential ages 
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of old passively evolving galaxies ( ISimon et al. 1 120051 ) . The second two-points subsample is 
determined from the high-quality spectra with the Keck-LR IS spectrograph of red-envelope 



galaxies in the redshift range 0.2 < z < 1 (jStern et al. Il2010[ ). The third subsample contains 



three data points obtained by using the 2-point correlation of SDSS lumino us red galaxies and 
takin g the BAO peak position as a standard ruler in the radial direction ( iGaztanaga et al. 
20091 ). These three points are considered as a direct measu rement of H(z) for t he first 
time and are independent of the earlier BAO measurements (jPercival et al. 1 12007| ) which 
constrains an integral of H(z) by using the spherically averaged (monopole) correlation. The 
parameters set can be determined by minimizing 



Xh(z) 



14 



[Hth{zi;p) - Hobsizi)f 



(9) 



When fully expanded, this expression includes Hq as a nuisance para meter. We margina lize 
over Hq usin g the Gaussian p rior Hq 
proposed in flWu fc Yu Il2007a[ l. 



73.8 ± 2.4 km s"^ Mpc'^ (IRiess et al. 11201 If ) as 



3. MODELS AND CONSTRAINT RESULTS 

For the sake of simplicity, we assume a spatially flat universe. We study several popular 
cosmological models. The chosen models, their parameters, and the abbreviations we use to 
refer to them are summarized in Tabled! In the following, we will outline the basic equations 
governing the background evolution of the universe in each model and give the observational 
results. 



Model 


Abbreviation 


Parameters 


Flat ACDM 


FA 




Flat constant w 


Fw 


Qrn, W 


Flat varying w (GPL) 


FGPL 


Qm, Wo, Wi 


Flat Dvali-Gabadadze-Porrati brane 


FDGP 




Flat Modified Polytropic Gardassian 


FMPG 




Flat Generalized Ghaplygin Gas 


FGGG 


As, 7 



Table 1: Summary of models. Note: since the parameter "a" has been used to describe the 
stretch-luminosity relation for the SNLS3, we utilize "7" to replace the original parameter 
"a" in the GGG model. 
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3.1. Flat ACDM 

The flat ACDM model is the simplest one that can explain the present accelerating 
cosmic expansion, and is generally considered to be the standard cosmological model. In 
this model, dark energy is a result of a nonzero cosmological constant. So, we have Qbe = 
Q\ = (1 — flm), with equation of state w = —1. The Friedmann equation in this case is 

= H^[QU^ + zf + (1 - QJ], (10) 

This model is consistent with almost all observations. Fig. [T] shows a comparison of different 
SN la data constraints on the model parameter. From this Figure, we see that irrespective 
of whether SALT2, SALT, or SIFTO light curve fitter is used, different SN la samples 
give fairly consistent results. However, this is not the case when the MLCS2k2 light-curve 
fitter is considered, different SN la sets lead to notably different results. For example, 
the best fit values are fi^ = 0.209lJ]:[!^^ ( ESSENCE), n^r,. = 0.324 l°-°i^g (Constitution) and 



= 0.403+^:^^;; (SDSS-II). in addition, as lSoUerman et al. I (120091), we also find that SDSS- 



II (SALT2) and SDSS-II (MLCS2k2) are clearly inconsistent. At the 95.4% confidence level, 
the constraints from SDSS-II (SALT2) and SDSS-II (MLCS2k2) have no overlap. In Fig. ^ 
we give the constraints from SNLS3 with other cosmological probes, from which we find 
that the main constraints come from the SN la+CMB/BAO and the H(z) data set has little 
effect. Using SNLS3+CMB/BA0+H(z), we obtain fi„ = 0.245l|5:[{^^ at the 68.3% confidence 
level. Compared with the value of Qm = 0.27 ± 0.04 obtained by Davis et al. (2007) from 
ESSENCE+BAO+CMB, we find that SNLS3+CMB/BA0+H(z) yields a tighter constraint 
on fi^ . 



3.2. Flat constant w 

This model is obtained by assuming that dark energy has a constant equation of state 
parameter w. Thus the Friedmann equation can be expressed as: 

= H^[Q^{1 + + (1 - QJ{1 + zf'+^\ (11) 

which depends on two free parameters, flm and w. Fig. [3] shows the constraints at the 68.3% 
confidence level from different SN la samples, from which we find that, for all SN la data, 
especially for the SDSS-II, different light-curve fitters still affect the results. In the SDSS-II 
case, there exists an obvious inconsistency between SALT2 and MLCS2k2 since there is no 
overlap at the 68.3% confidence level in this model. By fitting this model to the combined 
SNLS3 SN la, BAO/CMB and H(z) data, the 68.3% constraints are fi„ = 0.248l|]:JJ^^ and 
w = — I.O39I0J53. The contours of fi^ and w are plotted in Fig. HI w = — 1 is consistent 
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with observations at the 68.3% confidence level. As was the case for the ACDM, the H(z) 
has little effect on the derived cosmological parameters. 



3.3. Flat varying w (CPL) 



The degrees of freedom of the cosmic model increase when we allow the dark energy 
equation of state to vary with the cosmic time. In this case, the Friedmann equation for a 
varying w dark energy model is given by 

l+w(z') 



Hi 



+ zf + (1 - ^^„^)exp 



l + z' 



-dz' 



(12) 



In this paper, we consider the popular CPL ( IChevallier &: Polarski I l200ll : iLinder 1 12003[ ) 
parametrization w{z) = wq + Wiz/{1 + z). The above expression can then be simplified to 

?>WiZ 



Hi 



nm{i + zf + (1 - nm){i + ^)3(i+-o+-i) 



exp 



l + z. 



(13) 



This model has been investigated by using SNLS3 in Sll and lLi et al. I (120111 ). In Fig. 0, we 
show the marginalized 68.3% contours of Wq and Wi from different SN la data with different 
light-curve fitters. Consistent results are obtained and the tensions between different SN 
la sets and different light-curve fitters occurring in the ACDM and the constant w model 
disappear when the CPL is considered. However, this tension improvements are probably 
simply due to the fact that this model is not well constrained by current data sets, especially 
the Wi constraints are so bad. A combination of SNLS3+BA0/CMB+H(z) gives VLm = 
0.253l^;E^|, wq = -imrtoill, and wi = -0.344+^;^^^ at the 68.3% confidence level. The 
marginalized contours for Wi and Wq are plotted in Fig. [HI One can see that the ACDM [wq = 
— 1.0, Wi = 0.0) is included at the 68.3% confidence level for all observational data sets, which 
are consistent with what obtained in Sll from SNLS3+SDSS DR7 LRGs+WMAP7+ifo- 



3.4. Flat DGP model 



The DGP model (IDvali et al. Il2000l ). which accounts for the cosmic acceleration without 
dark energy, arises from a class of brane world theories in which gravity leaks out into the 
bulk at large distances. For a spatially fiat case, the Friedmann equation in the DGP model 
can be expressed as 

+ (14) 



H^ = Hl[^nUi + z] 

where f2m = 1 ^ 2 a/ ■ The parameter Tc represents the critical length scale beyond which 
gravity leaks out into the bulk and is related to this critical length by Qj-c = ^/(^'"c^o)- 
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The flat DGP model has the same number of free parameter as the ACDM. The constraints 
from different SN la samples are shown in Fig. [71 We obtain similar results to the ACDM 
for the DGP model, except that the latter favors a smaller Qrn- Fig- El shows the constraints 
on from SNLS3, CMB/BAO and H(z) data, from which we find that the SNLS3 gives 
a very small value for Qm and this value becomes large when the CMB/BAO data are 
added. The combined SNLS3+CMB/BA0+H('z) gives = 0. 2191°°^^, which is smaller 



than that obtained in IXu fc Wang I (120101 ): iLiang fc Zhu I (120101 ) from other observations, 
where = 0.2971°;°^^ and 0.285^°;^^^, respectively. 

It has been argued that the BAO data cannot be used to constrain the DGP model, since 
the details of structure formation are unclear in this model (iRydbeck et al. II2007I ). Moreover, 
it has been found that a tension between distance measures and horizon scale growth in 



the DGP exists a nd there is no way to alleviate it (ISeahra &: Hu I l2010l : ISong et al. 1 12007 



Fang et al. II2008I ). With these caveats in mind, we still present the CMB/BAO constraints 
for the sake of completeness. 



3.5. Flat Modified Polytropic Cardassian 



The Cardassian model (IFreese &: Lewis 1120021 ) is based on a modified Friedmann equa- 
tion in which an additional term Bp^ is added to the right hand side 



Pm + Bpl 



(15) 



where pm is the density of matter and n is a dimensionless parameter. This equation can be 
rewritten as 



8nG 



-Pr, 



( PCard \ 



\ Prr 



1-n 



-) 



(16) 



This model reduces to the fiat ACDM when = 0, and is related to the constant w model 
by n = 1 + w. Therefore, the results for the original Cardassian model can be directly 
transcribed from those of section 13.21 and there is no need of additional fit for that model. 
Here, we consider a modified Cardassian mod el, the modified po lytropic Cardassian model, 
in which the cosmic evolution is governed by (jWang et al. 1120031 ): 



-Pr^ 



f PCard \ 
^ Pm 



q{l-n) 



1/9 



(17) 



The above expression can be reexpressed as 

= h^{qu^ + zf[l + {Q;^' - i)(i + zf^(--'^f'^ 
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ACDM is recovered for q = 1 and n = 0. The marginalized 68.3% contours in the q — n 
plane for different SN samples are shown in Fig. |9l As was the case for the varying w CPL 
model, here we find that the differences between the constraints arising from different data 
and different light-curve fitters are negligible, which may also be due to the weak constraints 
on the model parameters from these observations. The constraint results from SNLS3 along 
with other data are shown in Fig. [TOl The ACDM is well consistent with all observations at 
the 68.3% confidence level. The combined analysis gives V lrr, = 0.248l^- ^ fl q = 1.098+^;^^^ 
and n = 0.014lnQ^fi, which is consistent with what given in IPavis et al. I ( 120071 ) (see Fig. 4) 



and in lWang fc Wu I fl2009f ) {Qr. 



271+° °^^ n 



0.824: 



-0.750 
-0.6225 



n 



-0.091 



-0.331\ 
-1.908 J 



3.6. Flat Generalized Chaplygin Gas 



The Chaplygin gas model ( iKamenshchik et al. 1 120011 ) unifies dark matter and dark 
energy by invoking a background fluid with an equation of state p cx p~'^, with 7 = 1 in the 
basic model. Here, we consider a generalization of this model (the Generaliz ed Chaplygin 



Gas, or GCG model) by allowing 7 to take arbitrary, but constant, values (IBento et al. 
2OO2I ). For the GCG, the Friedmann equation is 



= H',[n,{l + z)' + (1 - Q,) X {A, + (1 - A,){1 + ^)3(i+^))TT^] 



(19) 



where fib is the present dimensionless density parameter of baryonic matter which is related 
to the effective matter density parameter by rim. = flh + (1 ~ ^fe) (l — AsY^^'^'^'^^ and the 
WMAP7 observation gives nbh"^ = 0.02246 jKomatsu et al. Iboilj ). Here h is the Hubble 



constant in unit of 100 Km • s ^ ■ Mpc ^ and we marginalize over it by using the iRiess et al 



Jioiiil Ho prior. Ag is a model parameter which relates the pressure p and energy density 



p of the background fluid: Ag = p/p'^ 



7 



corresponds to the case of the cold dark 



matter plus a cosmological constant. The 68.3% contours on Ag and 7 from different SN 
data sets are shown in Fig. [TTl The results are similar to the constant w model. Inconsis- 
tency appears between SDSS-Il (SALT2) and SDSS-11 (MLCS2k2) and there is no overlap 
between them at the 68.3% confidence level. The constraints from SNLS3 along with other 
cosmological probes are shown in Fig. [121 By fitting this model to the combined SNLS3 SN 
la, BAO/CMB and H(z) data, the constraints are A, = 0.810lg:[j| and 7 = OmetoUt We 
find that the standard Chaplygin gas model (7 = 1) is ruled out by SNLS3-I-CMB/BA0 and 
SNLS3-hCMB/BA0-hH(z) data at the 95.4% confidence level, while the ACDM is consistent 



with t h em at the 68.3% co n fidence lev e l. Th is is consistent with the results of iDavis et al. 
J2007h : iLiang et al. I J201oh : ILu et al. I J2009h : IWu fc Yu I J2007al lbll 
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MODEL TESTING USING MODEL SELECTION STATISTICS 



In this section, we discuss the worth of models by applying model comparison statistics 
such a s \'^/do f (do f: deg rees of freedom), goodness of fit (Go F), Akaike informatio n cri- 
terion (lAIC: llAkaike 19741 ) and Bayesian information criterion f lBIC: llSchwarz 19781 ). The 
X^/dof describes how well the model fits a set of observations. The GoF simply gives the 
probability of obtaining, by chance, a data set that is a worse fit to the model than the 
actual data, assuming the model is correct. It is defined as 



GoF = T{u/2, x'm/T{u/2), 



(20) 



where F is the incomplete gamma function and u is the number of degrees of freedom. For 
a family of models, the best fit one has a minimum /dof value, while it has a maximum 
value of GoF. 

For a given data set, the candidate models may be ranked according to their AIC values, 
which can be calculated by 

AIC = -21n£ + 2A; (21) 

where C is the maximum likelihood, k is the number of model parameters. Given a set of 
candidate models for the data, the one which has the minimum AIC value can be considered 
the best. Then the relative strength of evidence for each model can be judged by using the 
differences (AAIC) between the AIC quantities of the rest of models and that of the best 
one. The models with < AAIC < 2 are considered substantially supported, those where 
4 < AAIC < 7 have less suppor t, while models with AAIC > 10 are essentially unsupported 
with respect to the best model (ISzydlowski &: Kurek II2008I ). 



The BIG, very similar to the AIC, is also a criterion for model selection among a finite 
set of models. It is defined as 

BIC = -2\nC + k\nN (22) 

where N is the number of data points used in the fit. The model which minimizes the BIG is 
the best fit one. As the AIC, the differences between the BIG of the rest of models and that 
of the best one (ABIC) is used for the judgement of the model, that is, < ABIC < 2 is 
considered as a weak, 2 < ABIC < 6 as a positive, 6 < ABIC < 10 as a strong and ABIC > 



10 as a very strong evidence favoring a better model (ILiddle I l2004l : ISzydlowski fc Kurek 



20081 ) . Both the AIC and BIG penalize the case of adding the model parameter to increase 
the likelihood through the introduction of a penalty term, which depends on the number of 
parameters in the model. But the coefficients in this penalty term are different for the AIC 
and BIG. It should be noticed that in the limit of large data (large N) the AIC tends to 
favor models with more parameters while the BIG tends to penalize them (jParkinson et al 
2005h . 
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The results for these model tests are shown in Table [2], [3l |H From these tables, we 
find that, the DGP model is favored by the SNLS3 data alone. However, if additional 
constraints such as CMB/BAO and H(z) are included, this is no longer true, and the DGP 
model is disfavored by most criteria. Instead, the fiat ACDM becomes preferred. The GCG 
and Fw are hardly distinguishable by most selection statistics and they become the favored 
non-standard models (excluding ACDM) when the CMB/BAO and H(z) are added. The 
variable w (GPL) and Gardassian (MPG) models are penalized the most because of their 
large number of cosmological parameters. 



5. Conclusion 

In this paper, we compare the constraints from the enlarged first three years data of the 
Supernova Legacy Survey (SNLS3) with those from other SN la samples. Several popular SN 
la samples, such as ESSENCE, Union2, Constitution and SDSS-II are used and the effects 
of different light-curve fitters on results are considered. We also discuss the observational 
constraints from the SNLS3, together with other two cosmological probes, CMB/BAO and 
Hubble parameter versus redshift. The ACDM and five nonstandard cosmological models are 
considered. We find that, for models with fewer free parameters (FA, FDGP, Fw and FGCG; 
see Tabled]), different SN la samples give fairly consistent results when analyzed with the 
SALT2, SALT, or SIFTO light-curve fitters, but this is not true when the MLCS2k2 light- 
curve fitter is used. Moreover, we find significant ten sion between SDSS-II s ample results 



when analyzed with SALT2 or MLCS2k2, as noted in ISoUerman et al. I (120091 ). The incon- 
sistencies between different samples and different light-curve fitters seem to be reduced for 
models with more free parameters (FCPL and FMPC). This improvement may be ascribed 
to the fact that current observational data has a weak constraint on the model parameters. 

By combining the SNLS3 with CMB /BAO and Hubble parameter versus redshift data, 
we find that, except for the fiat varying w model, the main constraints on model parameters 
come from the SNLS3-I-CMB/BA0 and the H(z) data has little effect. Moreover, we study 
the worth of models by applying model comparison statistics such as x^/dof, goodness of 
fit (GoF), Akaike information criterion and Bayesian information criterion. We find that 
the DGP model is the best one for SNLS3 alone while the fiat ACDM is preferred when 
the CMB/BAO or H(z) data is included. When only the nonstandard models are consid- 
ered (ACDM is excluded), the GCG and Fw are preferred by both SNLS3+CMB/BA0 
and SNLS3+CMB/BA0-|-H(z) data. These results are independent of the test methods, 
that is, different test methods give consisten t conclusions. Our res ults are slightly differ- 



ent from that of the analysis of the SDSS-II flSoUerman et al. Il2009[ ) and ESSENCE SN la 
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data flDavis et al. IboOTf ). In ISoUerman et al. I (120091 ). both SDSS-II (MLCS2k2) and SDSS- 
II (MLCS2k2) plus CMB/BAO prefer the DGP model, while, once t he SDSS-II (SALT2) 



is used, the ACDM is always the best fit one. In iDavis et al. I ( 120071 ). the combination of 



ESSENCE (MLCS2k2), BAG, and CMB favors the ACDM. Finally, the FCPL and FMPC 
suf fer in BIC test due to the extra model p arameters, which agrees with the conclusions 



in (IDavis et al. 1120071 : ISoUerman et al. Il2009f l. 
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Model X^/dof GoF(%) AAIC ABIC 

FDGP 419.464/472(0.8887) 96.04 0.00 0.00 

FA 419.691/472(0.8892) 95.98 0.23 0.23 

Fw 419.282/471(0.8902) 95.81 1.82 5.98 

FGCG.... 419.342/471(0.8903) 95.79 1.88 6.04 

FMPC.... 418.472/470(0.8904) 95.77 3.01 11.32 

FCPL 418.664/470(0.8908) 95.71 3.20 11.51 



Table 2: Summary of the model test results for the SNLS3. The models are listed in the 
increasing order of the AAIC values (please refer to Table [1] for the definitions of the model 
acronyms). When only the SNLS3 is considered, the fiat DGP is preferred. 



Model X^/dof GoF(%) AAIC ABIC 

FA 422.741/474(0.8919) 95.61 0.00 0.00 

FGCG.... 422.487/473(0.8932) 95.38 1.75 5.91 

Yw 422.603/473(0.8935) 95.34 1.86 6.02 

FCPL 421.124/472(0.8922) 95.52 2.38 10.71 

FMPC... 421.213/472(0.8924) 95.49 2.47 10.79 

FDGP 429.995/474(0.9072) 92.71 7.25 7.25 



Table 3: Summary of the model test results for SNLS3+CMB/BA0. The models are listed 
in the increasing order of the AAIC values (please refer to Table [1] for the definitions of the 
model acronyms). When the additional CMB/BAO is included, the fiat DGP is disfavored 
and the fiat ACDM becomes preferred. 
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Fig. 1. — The constraints on the fiat ACDM model from several popular SN la data sets 
compiled with different light-curve fitters. "Const" represents for the Constitution SN la 
sample. 




Fig. 2. — The constraints on the flat ACDM model. The blue dashed, red dotted and solid 
lines represent the results from SNLS3, SNLS3+CMB/BA0 and SNL3+CMB/BA0+H(z), 
respectively. 
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Fig. 3. — The 68.3% contours for the flat constant w model from several popular SN la data 
sets compiled with different light-curve fitters. "Const" represents for the Constitution SN 
la sample. 
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Fig. 4. — The 68.3% and 95.4% confidence level contours for the flat constant w model. Left 
panel shows the constraints from SNLS3, CMB/BAO and their combination, and right panel 
is the case with the H(z) data included. The red star {0.^ = 0.25, w = —1.0) represents the 
ACDM with the best fit Qrn 
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Fig. 5. — The marginalized 68.3% contours of Wq and Wi in the flat CPL model from 
several popular SN la data sets which are compiled with different light-curve fitters. "Const" 
represents for the Constitution SN la sample. 
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Fig. 6. — The marginalized 68.3% and 95.4% confidence level contours on wo — wi for the 
varying w model with the CPL parametrization. Left panel shows the constraints from 
SNLS3, CMB/BAO and their combination, and right panel is the case with the H(z) data 
included. The red star {wq = —1.0, Wi = 0.0) represents the ACDM. 
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Fig. 7. — The constraints on the flat DGP model from several popular SN la data sets 
compiled with different light-curve fitters. "Const" represents for the Constitution SN la 
sample. 




Fig. 8. — The constraints on the flat DGP model. The blue dashed, red dotted and solid 
lines represent the resuhs from SNLS3, SNLS3+CMB/BA0 and SNL3+CMB/BA0+H(z), 
respectively. 
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Fig. 9. — The marginalized 68.3% contours for the flat modified polytropic Cardassian model 
from several popular SN la data sets which are compiled with different light-curve fitters. 
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Fig. 10. — The marginalized 68.3% and 95.4% confidence level contours for the flat modified 
polytropic Cardassian model. Left panel shows the constraints from SNLS3, CMB/BAO 
and their combination, and right panel is the case of H(z) data included. The vertical line 
(g = 1.0) represents the constant w{= n — 1) model, and the red star (g = 1.0, n = 0.0) 
corresponds to the ACDM one. 
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Fig. 11. — The constraints on the flat GCG model from several popular SN la data sets which 
are compiled with different light-curve fitters. "Const" represents for the Constitution SN 
la sample. 
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Fig. 12. — The 68.3% and 95.4% confidence level contours for the fiat GCG model. Left 
panel shows the constraints from SNLS3, CMB/BAO and their combination, and right panel 
is the case with the H(z) data included. The horizontal line (7 = 0.0) represents the stan- 
dard cosmological constant model (ACDM) and 7 = 1.0 corresponds to the original Chap- 
lygin gas one. The original Chaplygin gas model is ruled out by SNLS3+CMB/BA0 and 
SNLS3+CMB/BA0+Hubble at the 95.4% confidence level, while, except for the CMB/BAO 
data, the LCDM is allowed by observations at the 68.3% confidence level. 
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Model X^/dof GoF(%) AAIC ABIC 

FA 436.106/488(0.8937) 95.57 0.00 0.00 

FGCG.... 435.931/487(0.8951) 95.31 1.83 6.02 

Yw 435.938/487(0.8951) 95.31 1.83 6.02 

FMPC... 435.656/486(0.8964) 95.08 3.55 11.93 

FCPL 435.900/486(0.8969) 94.99 3.79 12.18 

FDGP 443.030/488(0.9078) 92.85 6.92 6.92 



Table 4: Summary of the model test results for SNLS3+CMB/BA0+H(z). The models are 
listed in the increasing order of the AAIC values (please refer to Tabled] for the definitions of 
the model acronyms). When the additional H(z) is included, the flat DGP is still disfavored 
and the flat ACDM remains preferred. 



